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Abstract 

CO 
O . 

We develop a spectral representation formalism to calculate the Casimir force, in the 
Qh non-retarded limit, or van der Waals force between a spherical particle and a substrate, 

both with arbitrary local dielectric properties. The spectral formalism allows to study 
the system as a function of its geometrical properties separately from its dielectric 
properties. The calculated force is attractive, and at a small separations it is orders of 



magnitude larger for nanometric-size spheres than for micrometer particles. We also 
found that the force depends more on the dielectric properties of the sphere than of 
the substrate. 



Recent advances in micro and nano devices have opened the possibility of studying quantum 
phenomena that occur at these length scales. Such is the case of the Casimir force P that is 
a macroscopic manifestation of the quantum vacuum fluctuations, as predicted by quantum 
electrodynamics. The textbook example [21 EH E] consists of two parallel neutral conducting 
plates which attract each other. The first experimental measurements were done in 1951 
using dielectric materials [S], and in 1958 using conductors jBJ. These measurements have 
large errors, and up to recently, it was possible to perform measurements with about 15% of 
precision on truly parallel metal surfaces [Zj . The difficulty of keeping the two plates parallel 
at separations of few nanometers makes it easier to measure the Casimir force between a 
sphere and a plane [HI El QUI EH E2] • In this case, the Casimir theory for parallel plates 
can be extended using the proximity theorem j^j. The approximation is valid when the 
minimum separation between the sphere and the plane is much smaller than the radius of 
the sphere. This theorem was employed to corroborate experimental measurements of the 
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Casimir force between a plane and a large sphere [TU1 ITT| IT2] . However, it is well known 
that quantum effects become more evident as the size of the system decreases. Thus, the 
question of how important are the Casimir effects on nanometric-size spheres is still an open 
question of fundamental importance. In 1948 Casimir and Polder ^3] calculated the force of 
a polarizable atom near a perfect conductor plane considering the influence of retardation, 
and finding a correction to the London or van der Waals forces. Retardation effects are 
important if we consider that the distance between the atom and the plane is larger than 
the characteristic length of the system. Complementary theories are necessary to handle 
nanometer-size systems with real dielectric properties. Within this context Ford cal- 
culated the force between a perfectly conducting wall and a sphere with a Drude dielectric 
function. After a delicate cancelation of terms in the equations, he obtained a force that 
changes from attractive to repulsive in an oscillatory fashion depending on the relative dis- 
tance between the sphere and the surface. However, this oscillatory behavior has not been 
observed experimentally [111311011111112], and has not been predicted by other theories [5^IT5]. 

In this work, we develop a spectral representation formalism to calculate the force between a 
sphere and a substrate. The advantage of this spectral representation is that we can separate 
the contribution of the dielectric properties of the sphere and substrate from the contribu- 
tion of its geometrical properties. Since results for large spheres [5] and large distances [T3] 
are known, we restrict ourselves to the case of nanometric-size spheres and distances of few 
nanometers. In this case, it is not necessary to consider retardation effects, therefore, we 
work in the quasi-static limit such that the radius of the sphere and the minimum sepa- 
ration between the sphere and the plane, are smaller than the characteristic length of the 
system [15 . In this regime, the Casimir force is commonly known as the van der Waals or 
London force |3]. 
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Figure 1: Schematic model of the system. 

We consider a homogeneous sphere of radius R, electrically neutral and with a local dielectric 
function e sp h(^). The sphere is suspended at a minimum distance z above a substrate (see 
Fig. 1) which is also neutral and has a local dielectric function e su b(u;). The space or ambient 
between the sphere and substrate is vacuum (e am b = 1). The quantum fluctuations of the 



electromagnetic field induce a polarization in the sphere which can be described by a point 
dipole located at its center, 

Psph H = a(^)E vac H, (1) 

where a(u) = [e sp h(co>) — e am b]/[e sp h(c<j) + 2e am b]-R 3 , is the polarizability of the sphere that is 
assumed to be polarized uniformly jTHj, and E vac (u) is the electromagnetic field associated 
to the vacuum fluctuations. We can rewrite the polarizability as 

n — u{(jj) 

where n = 1/3 is a constant and u{uj) = [1 — e sph (co')/e amb ]~ 1 is a variable that only depends 
on the dielectric properties of the sphere and the ambient. When the sphere is near a 
substrate, it induces a charge distribution on the substrate that can be seen as a dipole 
image, such that 

PsubM = /c(w)M • PsphO), (3) 

where the term / C M satisfies the boundary conditions of the system, being M a diagonal 
matrix whose elements depend on the choice of the coordinate system, and f c {uj) = [e amb — 
e S ub(^)]/[eamb + e su b(^)] a the contrast factor that only depends on the dielectric properties 
of the substrate and the ambient. The induced charge distribution on the substrate produces 
a field which also modifies the sphere's dipole moment through a local field. Thus, the total 
induced dipole-moment on the sphere is 

p sph M = a(u) [E va » + T • p suh (oo)} , (4) 

where T is the dipole-dipole interaction tensor and, in the non-retarded limit, it takes the 
form 

m 3f r — 1 

T = — ^> (5) 

where 1 is a unitary matrix, r = (0, 0, 2(z + R)) is the vector from the center of the image 
dipole to the center of the sphere, r = r/r, and r = |r|. Given the symmetry of the system the 
diagonal components of M are (—1, — 1, 1), and there are only three independent components 
of T, one perpendicular to the surface plane and two parallel to this plane. 

The frequencies that satisfy the boundary conditions of the system are those at which the 
sphere is polarized. These frequencies are known as the proper electromagnetic modes of the 
system, and we denote them like u> s . Then, the total energy of the system is £ = J2 S l/2hu s . 
A convinient way of determining these proper electromagnetic modes is using a spectral 
representation formalism that we derive as follows. 

First we rewrite Eq. (@J) using the expression of the polarizability from Eq. (J2J), as 

[-u( W )l + H]. Psph ( W ) = V™ (6) 

where V vac = n -R 3 E vac , and H = n [l — f c (uj)R 3 T'M] is a dimensionless matrix that only 
depends on the geometry of the system. To find the solution of Eq. ©, consider the case 
when f c (u>) is real, then EI is a real and symmetric matrix. In this case, we can always find 



a unitary transformation that diagonalizes it, U 1 HU = n s , being n s the eigenvalues of EL 
Furthermore, the solution of Eq. JBJ) is given by 

Psph M = G( M )V vac , (7) 

where G(u) = [— u(u)l + EI] -1 is a Green's operator. The r/'th element of G(u) can be 
written in terms of the unitary matrix U as ^7j 



u — n s 

s 

The poles of G(u), that is u{uj) = n s , give the frequencies of the proper electromagnetic 
modes, lo s , of the system ^H]- We now calculate the Casimir interaction energy as the 
difference between the energy when the sphere is at a distance z from the substrate and the 
energy when z — > oo, that is, 

^Ef^T- ( 9 ) 

s s' 

The eigenfrequencies uj s i are obtained from the poles of Eq. (jHJ) when z — > oo, or by substi- 
tuting e su b(co | ) = e am b in f c . Note that it is not necessary to do any renormalization or any 
delicate cancelation to calculate the energy. Alternatively, we can also find the density of 
states using the Green's function definition and then calculate the energy of the system, as 
we show in the appendix. 

The advantage of the spectral representation is that we can separate the contribution of the 
dielectric properties of the sphere from the contribution of its geometrical properties. As we 
mentioned, the material properties of the sphere are contained in the spectral variable u, 
while the geometrical properties of the system, like the radius of the sphere and the separation 
of the sphere to the substrate are in the matrix HI. Furthermore, EI is a dimensionless matrix 
that depends on the ratio z/R. Its eigenvalues are independent of V vac and of the dielectric 
properties of the sphere. And the dielectric properties of the substrate are in f c which is a real 
function even for dispersive materials A similar spectral representation was proposed 
years ago to study the effective dielectric properties of granular composites [20J. The results 
discuss here are calculated as follows. First, we construct the matrix EI for a given z/R, and 
we diagonalize it numerically to find its eigenvalues n s . Considering an explicit expression 
for the dielectric function of the sphere, we calculate the proper electromagnetic modes uj s 
trough the relation u(uj s ) = n s . Once we have u s , we calculate the energy according with 
Eq. ©. Here, we use the Drude model, such that e sp h(co>) = 1 — ujp/[uj(uj + z/r)], where uj p 
is the plasma frequency and r is the relaxation time. We present results for potassium (K), 
gold (Au), silver (Ag) and aluminum (Al) spheres with fkv p = 3.80, 8.55, 9.60, and 15.80 eV, 
and {rujp)^ 1 = 0.105, 0.0126, 0.00188, and 0.04, respectively. We have considered substrates 
whose dielectric function is real and constant in a wide range of the electromagnetic spectrum 
as sapphire (A1 3 2 ), and titanium dioxide (Ti0 2 ), with e su b = 3.13, and 7.81, respectively. 
Then, the corresponding contrast factors are f c = -0.516, and -0.773. We have also considered 
the case of a perfect conductor substrate (denoted by Inf) with e su b — > oo and f c — — 1. 



In Fig. 2, we show the energy as a function of z/R. In general, we observe that the energy 
shows a power law of (z/R)~ 3 . This behavior is independent of the material properties, 




Figure 2: Energy as a function of z/R. Each panel shows the results for different substrates 
(perfect conductor, AI3O2, and TiC^). 



and it is inherent to the dipole-dipole interaction model. This is consistent with the result 
found by Casimir and Polder [IHJ for a polarizable atom, and with the measurements by 
Mohideen et al. but it is contrary to the oscillatory behavior calculated by Ford |14j . 
The value of the energy varies with the substrate, for example, at small distances it is about 
two times larger for a perfect conductor substrate than for AI3O2, while the TiC>2 case is 
between them. This is easily explained if we look at the contrast factor values for each 
substrate, where one can see that as f c — > —1, the energy is larger. For all the substrates, we 
found that V becomes larger as the plasma frequency of the metal also does. In conclusion, 
we found that V is large when f c — > — 1 and u p is large, recovering the limit for perfect 
conductors. Therefore, the energy is largest (smallest) for an Al (K) particle over a perfect 
conductor (A1 3 2 ) substrate. When the sphere is at a distance larger than 2R, the energy 
is very similar, independently of the dielectric properties of the sphere and substrate. 

Let us first analyze the force as a function of the geometrical properties, that is, as a function 
of R and z. In Fig. 3 we show the Casimir force calculated as 



F = — - 



d£ 

dz 



1 d£ 
Rdjz/Rj' 



(10) 



In all cases, we obtain an attractive force such that, as R is smaller the force increases. When 
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Figure 3: Casimir force as a function of z. In the left (right) panel we show results for a 
K (Al) sphere over an A1 3 2 (Inf) substrate. The different curves correspond to spheres of 
10 nm, 50 nm, 100 nm and 500 nm of radius. 

the sphere is almost touching the substrate (z ~ nm) the force is ten times larger for a 
sphere of R = 10 nm than the one of R = 100 nm, and increases fifty times for a sphere 
of R = 500 nm. As a function of z, the force for the sphere with R = 500 nm seems to 
be almost constant from to 40 nm compared with the other curves. This is an artifact of 
the scale since all curves have a power law of z~ 4 , and they are proportional to R 2 . This 
implies that for a distance z < 10 nm the force is larger for the smallest sphere; however, 
at a larger distance the force is larger for larger spheres, while for very large distances, the 
force is independent of R. Furthermore, the force for the sphere with R = 10 nm decreases 
about three orders of magnitude as the separation of the sphere goes from to 40 nm, 
independently of the dielectric properties of the system. On the other hand, with the proper 
combination of dielectric functions of the sphere and substrate it is possible to modulate the 
magnitude of the Casimir force. Here, we show the force for an Al sphere over a perfect 
conductor which is one order of magnitude larger than the force between the K sphere over 
A1 3 2 , as it is expected. 

From Fig. 4 we analyze the force as a function of the dielectric properties of the particles 
and substrate. In all cases, we found the same dependence of the force as a function of 
z, independently of the dielectric functions of both, sphere and substrate. Although the 
dielectric function of the substrate is important, the dependence on the dielectric function of 
the sphere is more critical in the magnitude of the Casimir force. Indeed, the force is larger 
for increasing values of uj v . In particular, we found that the force for an Al sphere is almost 
ten times larger than for a K sphere. On the other hand, for a given sphere the Casimir 
force increases at most by a factor of three, when the substrate is changed from A1 3 2 to a 
perfect conductor. 



The spectral representation formalism allows us to calculate the force between a sphere and a 
substrate in a range of sizes and separations where the proximity theorem is not applicable. 
For example, for a Au sphere of radius 100 nm on top of perfect conducting plane, the 
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Figure 4: Casimir force as a function of z for spheres of K, Au, Ag, and Al with R = 50 nm 
over AI3O2, and Ti02- 

value of the Casimir force using the proximity theorem yields a force about three orders of 
magnitude of the values obtained in this paper. This is due to the linear dependence of 
the proximity theorem with the radius of the sphere and that the geometrical effects of the 
sphere are not included. On the other hand, systems where the proximity theorem is used, 
such as the experiments by the group of Mohideen and collaborators [TTJ H2] do not employ 
homogeneous metallic spheres, rather coated dielectric spheres, making it difficult to employ 
the spectral representation formalism. 

In conclusion, we have developed a spectral representation formalism within the van der 
Waals approximation to calculate the Casimir force between a sphere and a substrate. This 
spectral formalism separates the geometrical properties contributions from dielectric prop- 
erties contributions on the Casimir effect in the non- retarded limit. We found that at very 
small distances, the force can increase orders of magnitude as the size of the particle becomes 
smaller. We have also observed that the correct choice of the dielectric properties of both, 
sphere and substrate, can increase or decrease the force by orders of magnitude. 

This work has been partly financed by CONACyT grant No. 36651-E and by DGAPA-UNAM 
grants No. IN104201 and IN107500. 



A Appendix 

By using the Green's function from Eq. (8), one can find the density of states of the system 
as a function of the spectral variable u, as p(u) = — -Im [TrG(w)] . In the particular case, 
when the sphere is metallic and its dielectric function is described by the Drude model, 
u{uj) = uj{u + ir)/ujp. Then, the density of states as a function of the frequency takes the 



following explicit form 

7T ^ \(CJ 2 - UJplsY + {uJ/tYJ 

and the energy can be calculated as V = \ J °° ftw p(uj)duj . 
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